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Abstract 

In a recent paper by Hooper and Goodenough, data from the Fermi Gamma Ray Telescope was 
analyzed and an excess of gamma rays was claimed to be found in the emission spectrum from 
the Galactic Center Region. Hooper and Goodenough suggest that the claimed excess can be 
well explained by 7-10 GeV annihilating dark matter with a power law density profile if the dark 
matter annihilates predominantly to tau pairs. In this paper we present such a dark matter model 
by extending the MSSM to include four Higgs doublets and one scalar singlet. A 7L<i symmetry is 
imposed that enforces a Yukawa structure so that the up quarks, down quarks, and leptons each 
receive mass from a distinct doublet. This leads to an enhanced coupling of scalars to leptons and 
allows the model to naturally achieve the required phenomenology in order to explain the gamma 
ray excess. Our model yields the correct dark matter thermal relic density and avoids collider 
bounds from measurements of the Z width as well as direct production at LEP. 
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I. INTRODUCTION 



Recently, Hooper and Goodenough examined the first two years of Fermi Gamma Ray 
Space Telescope (FGST) data from the inner 10° around the Galactic Center They found 
that the gamma ray emissions coming from between 1.25° and 10° of the Galactic Center 
is consistent with what is expected from known emission mechanisms such as cosmic rays 
colliding with gas to produce subsequently decaying pions, inverse Compton scattering of 
cosmic ray electrons, and known gamma ray point sources. In order to model the gamma ray 
background within 2° of the Galactic Center, Hooper and Goodenough model the emission 
of the Galactic black hole Sgr A* as a power-law extrapolated from higher energy HESS 
observations. Comparing the FGST measurements to this background, Hooper and Good- 
enough found that it agrees very well with FGST data between 1.25° — 2° but found an 
excess in the observed gamma ray intensity within 1.25°. It has been pointed out by Ref. 

however, that a simple power-law extrapolation of HESS data may understate the flux of 
the central point source Sgr A* as the slope of its spectrum may deviate from the constant 
HESS results below an energy of 100 GeV. 

The authors of Ref. |l| showed that the increased gamma ray emissions are well described 
by annihilating dark matter that has a cusped halo profile (p oc r -7 , with 7 = 1.18 to 1.33) 
provided that the dark matter satisfies three basic conditions. The conditions required of the 
dark matter are 1) that it have a mass between 7 — 10 GeV, 2) that it annihilate into r-pairs 
most of the time, but into hadronic channels 15 — 40% of the time, and 3) that its total 
annihilation cross section yield a thermal average within the range (av) = 4.6 x 10~ 27 — 
5.3 x 10 -26 cm 3 /s. It should be noted that the results of Hooper and Goodenough are 
controversial, and the Fermi-LAT collaboration itself has not yet published official results. 
In addition, other background related explanations for the gamma ray excess have been 
proposed such as the existence of a pulsar near the Galactic Center [3|. In this paper we 
proceed with the assumption that the analysis of Hooper and Goodenough is correct. The 
astrophysical and particle physics implications of this finding are discussed in Refs. 0, 0]. 

In this paper we construct a dark matter model satisfying the above conditions by adding 
a singlet to the supersymmetric leptophilic Higgs model (SLHM) jg]. In the SLHM the up 
quarks, down quarks, and leptons, each receive mass from a separate Higgs doublet. For 
our purposes, the salient characteristic of the SLHM is that it endows the leptons with an 
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enhanced coupling to one of the scalars. This provides a natural mechanism for dark matter 
particles to annihilate predominantly into r-pairs. This model of dark matter is able to 
successfully account for the FGST observations, yields the correct relic density, and evades 
relevant collider bounds such as measurements of the Z width and direct production at LEP. 
The idea of a leptophilic Higgs has been studied as a possible explanation for the e ± excess 
observed by PAMELA and ATIC in Ref. \. However, this entails a 100 GeV - 1 TeV dark 
matter particle, while our model requires a light, O(10) GeV dark matter particle. There 
also exist some other models that can explain the Galactic Center gamma ray excess s|. 

In addition to explaining the FGST observations, such a model of light dark matter is 
also capable of describing observations by the CoGeNT [9] and DAM A collaborations 10 ]. 
CoGeNT has recently reported direct detection signals that hint at the presence of O(10) 
GeV dark matter compatible with the light dark matter interpretation of DAMA's annual 
event rate modulation. Ref. ll| showed that dark matter with a mass between 7 — 8 GeV 
that has a spin independent cross section approximately between a si = 1 x 10~ 40 — 3 x 
10~ 40 cm 2 is consistent with both CoGeNT and DAMA signals. Although the XENON [12 1 
and CDMS jl3| collaborations challenge this report, Ref. {5] has pointed out that "zero- 
charge" background events lie in the signal region. The authors suggest that the bound 
could possibly be loosened if a modest uncertainty or systematic error is introduced in the 
energy scale calibration near the energy threshold. Although our model is able to explain 
the reported observations of the CoGeNT and DAMA collaborations, it is not dependent 
upon their validity. By simply moving to another region of parameter space our model can 
coexist with the absolute refutation of CoGeNT and DAMA while continuing to explain the 
FGST results and avoiding collider bounds. 

Our paper is organized as follows. In Section [Til we introduce the setup of the model and 
calculate the mass matrices for the scalars and the neutralinos. In Section [TTT1 we describe the 
process by which the dark matter annihilates into Standard model particles and calculate 
the relevant cross sections for a benchmark point in parameter space. We also show that 
the resultant relic density is consistent with current cosmological measurements. In Section 
IIVI we discuss possible direct detection and in Section [V] we discuss relevant bounds for 
this model and show that it is currently viable. Lastly, we conclude with Section I VI I and 
summarize the results of the paper. 
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II. THE MODEL 



In this model the quark and lepton content is that of the MSSM. To this we add four 
Higgs doublets, H u , H d , H , and He, with weak hypercharge assignment +1/2, —1/2, +1/2, 
and —1/2 respectively. The third Higgs doublet is necessary to achieve a leptonic structure, 
while the fourth doublet is required for anomaly cancelation. In order to avoid problems 
with the Z decay width, we introduce a singlet S that acts as O(10) GeV dark matter. The 
idea of adding a light singlet to the MSSM to act as dark matter was also considered in 
[jjj], while the use of a singlet for other purposes such as solving the /i problem was first 

is given by 

W = yJJQEu - y d DQH d - y e ELHe + ii q H u H d + ^H H e 

^1^1^ (2-1) 

+ K q SH u H d + K e SH Ht + XfS + - X 2 S 2 + - k s S 3 , 

where the hats denote superfields. In the superpotential we introduced a Z 2 symmetry 
under which if , He and E are odd while all other fields are even. The symmetry enforces 
a Yukawa structure in which H u gives mass to up-type quarks, H d to down-type quarks, 
and He to leptons, while Hq does not couple to the quarks or leptons and is called the inert 
doublet. It is introduced to ensure anomaly cancellation. The Z2 symmetry is broken in 
V^oft so that we have: 1 

V so{t = ml\H u \ 2 + m 2 d \H d \ 2 + m 2 \H \ 2 + m]\He\ 2 + m 2 s \S\ 2 
n\H u H d + nlH Hi + ii\H u He + /J.lH H d 



(2.2) 

+ ^ a SH u H d + fibSH He + fi c SH u He + fi d SH H d 

+ m 2 u0 HlH + m 2 dl H\He + t 3 S + b 2 s S 2 + a s S 3 + h.c.) . 

The breaking of the Z2 symmetry is discussed in greater detail in Appendix |A] The Higgs 
sector potential is given by V — Vd + Vf + Koft- Letting a a denote the Pauli matrices for 
a = 1,2,3, the D-term is simply 

,2 , 2 



V D = 9 - \Hl<y a H u + H\a a H d + # V a # + H\a a H, 



a 



a' 2 1 

_j_ y I tt 1 2 I tt 1 2 I I tt |2 I rr |2 

+ -3- M-nul — l-nrfl + l-nol — |^| 
o I 



(2.3) 



1 In Ref. [6( the soft breaking terms m1 Q HlH + m de H^He + h.c. were omitted. 



where g and g' are the SU(2) and U(l) gauge couplings respectively. The F-term and V so ft 
combine with the D-term to yield the following potential 



V 



(fi 2 q + m 2 u )\H u \ 2 + (fj, 2 q + ml)\H d \ 2 + (rf + m 2 )\H \ 2 + + mj)\H t \ : 
(fj,l + K g \l)H u H d + (i4 + KiX\)H H t + f/ 3 H u Ht + nlH H d + h.c. 



+ 



+ 



- (m 2 u0 HlH + m 2 de H d H e + h.c.) + (m 2 s + X 2 2 ) \S\ 2 
(t 3 + X 2 1 X 2 )S+ (b 2 s + k s X 2 2 )S 2 + a s S 3 + h.c] + k s X 2 \S\ 2 (S + S*) + k 2 s \S\ 
fi a (H u H d ) S + fi b (H Hz) S + fi c (H u H e ) S + fi d (H H d ) S + h.c." 



(2.4) 



+ A 



K q {H u H d ) +k 1 {H H 1 ) 



S* + K s 



,(H u H d ) + K t (H H t ) (S 2 )* + h.c 



{ Kq fi q (\H u \ 2 + \H d fj + W (|# | 2 + TO 2 ) } {S + 5*) 



+ K\ \E U Y + |^| 2 ) \S\ 2 + k] ( |tf | 2 + \Htf )\S\* + V D . 



The singlet S acquires the vev (S) = v s /\/2 while the Higgs doublets acquire the vevs: 



TO) 




TO 



v/2 




.to 




TO 




(2.5) 



Letting v 2 w = v 2 u + v 2 d + v 2 + vj so that v 2 cw = AM 2 z /(g 2 + g' 2 ) « (246 GeV) 2 , we define the 
mixing angles a, /3, and fa by the relations tan/3 = v u /v d , tan fa = v /v£, and tan 2 a = 
(t> 2 + v d )/(v 2 + w 2 ). These definitions lead to the following parameterization of the Higgs 
vevs: 



(2.6) 



v u — v ew sin a sin fa w rf = t; cw sin a cos /3, 
= ^cw cos a sin = v cw cos a cos fa. 

In order to avoid increasing the Z width or violating other known bounds, we want 



the light dark matter to separate from the other neutralinos and be mostly singlino s, the 
fermionic component of the singlet S. This is accomplished by taking the parameters n q and 
Ki to be small, which eliminates most of the mixing between the singlino and the Higgsinos 
[see Eq. ( I2.10p ], It can then be easily arranged to have the singlino be the lightest of the 
neutralinos. A possible mechanism for explaining the small size of n q and k,£ is discussed 
in Appendix |A] Small values of n q and Kg also leads to reduced mixing between the scalar 
singlet and the Higgs doublets as can be seen from Eq. (I2.4p . A small amount of mixing 
is of course required since we desire the lightest scalar, which is mostly singlet, to couple 



to r-pairs in order for the dark matter to annihilate to r + r and other Standard Model 
particles. This mixing is generated by the soft supersymmetry-breaking parameters /i a , fib, 
He, and fi d . 

It is sufficient for n q and to be (9(10~ 2 ), which is what we use in our numerical 
calculations (see Table HI and HT|) . Though the scalar mass matrices are quite complicated 
in general, they simplify considerably in the limit of vanishing K q and Kg. The numerical 
calculations in the sections that follow have been determined using the general matrices, 
but for compactness we present only the simplified matrices here. In the {h u , hd, ho, hi, h s } 
basis, the neutral scalar mass matrix is given by 



M 2 N 



M 2 7ti 2 
lP T Ml 



(2.7) 



ss 



where the matrix M 2 is given by M 
are given by 



- M| LHM + AM 2 + AMf and the terms rrt 2 and M : 



ss 



and 



M 



ss 



3(a s + k s X 2 )v 2 s + 2\/2k 2 s vI - 2t 3 - 2A 2 A 2 + (n a v u Vd + HbVQVi + [i c v u vi + fidVoV d ) 

V2~v s 



The matrix M| LHM is the neutral scalar mass matrix from the ordinary SLHM, which can 
be found in while the matrices AM 2 and AM 2 are given by 



AM 2 



( 
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-V s fld 







-Vgfld 
-V s flb 



-V s /J> c 



j£ (H b v + fi c v u ) J 
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The pseudoscalar mass matrix, in the {a u , ad, do, ag, a s } basis, is similarly given by 

M 2 -Wt 2 
-7?t 2T M 2 



M 2 



(2.8) 



L ss 



where M 2 = M| LHM + AM 2 + AMf. The matrix M| LHM is the pseudoscalar mass matrix 
from the ordinary SLHM while AMf is the matrix obtained from AMf by changing the sign 
of every off-diagonal entry. Lastly, Mj s is given by 



Mis 



y/2'i 



fJ> a V u Vd + VbVQVi + [l c V u Vi + HdV Vd - 2A 2 A 2 



2t 3 - (9a, + k s X 2 )v 2 s - 4V2 (o 2 + k s X 



2! u s 



The chargino mass matrix, on the other hand, is rather simple even with nonvanishing K q 
and Kg. Letting h u , hj, ho, and he denote the Higgsino gauge eigenstates, the chargino mass 
matrix, in the {W + , h^,fiQ, W~ , h^,hj} basis, is given by 

\ 



M x± 



( 








M 2 


gvd 


gvt 











gv u 


» q + %v s 














gv 





h + % 


M 2 


gv u 


gvo 











gvd 





















H + f 2 v s 












(2.9) 



/ 

Like the chargino mass matrix, the neutralino mass matrix is simple. The neutralino mass 
matrix, in the {B°, W°, h u , hd, h , he, 5} basis, is given by 



M v 
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M x 





\g'v u 


-\g'v d 


\g'vQ 


-\g've 





\ 







M 2 


-\gv u 


\gvd 


-\gvo 


\gve 









\ g'v u 


-\gv u 













T2 Vd 






-\g'vd 


\gvd 



















\g' v o 


-\gv 











1*1 + %v. 








-\g've 


\gvi 

















V 










V2 Uu 




A2 


-f V2k s v s 


J 



(2.10) 

When n q and Kg are small, the singlino part of the above matrix separates from the wino, 
bino, and higgsinos, and the singlino mass can be well approximated by 

m xi « A 2 + V2k s v s . (2-11) 
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TABLE I: Benchmark Point A 
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TABLE II: Benchmark Point B 

The O(10) GeV LSP can be arranged with some tuning of the parameters in order to achieve 
a cancelation between A 2 and the product k s v s in Eq. (12. lip . Though the smallness of K q 
and Kg is technically unnatural, we remind the reader that a possible mechanism to make 
them small is discussed in Appendix |A] 

In the following sections, we calculate the relevant cross sections and quantities of interest 
using benchmark points A and B, found in Tables [I] and [TT] respectively. While both of 
these benchmark points can explain the Galactic Central region gamma ray excess, the spin 
independent direct detection cross section corresponding to benchmark point A lies within 
the region favored by CoGeNT and DAMA. In contrast, we will show that benchmark point 
B satisfies CDMS bounds that exclude CoGeNT and DAMA. Relevant quantities have been 
calculated for several additional benchmark points as well, and their values are summarized 
in Table ED of Appendix [Bl 
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FIG. 1: The dominant diagram of dark matter annihilation into fermions. Here a\ is the lightest 
pseudoscalar. 

III. ANNIHILATION TO FERMIONS 



In this section, we will show that this model can achieve the conditions needed to explain 
the gamma ray excess in the Galactic Center region. In order to calculate the dark matter 
cross section, we need the interactions between Higgs and fermions: 

C D 4= [h s ss — m s §7 5 s] 

v2 



K q 
2^2 



h u sh d - ia u sj 5 h d + h d sh u - ia d srf'h u + h.c. 



'h shi — iaos^he + h e sh — ia£S"y 5 h + h.c. 



2V2 



(3.1) 



f={u,d,e} j J 



(hffjfj ~ ia ffj7 5 fj) , 



where is the mass of the fermion fj, Vf is the vev of /-type scalars, and j runs over the 
fermion generations. In the limit n q , k,£ — > 0, the higgs-higgsino-singlino interactions vanish. 
We can expand (av) in powers of the dark matter velocity squared v 2 : 

( av ) = a + bv 2 + .... (3.2) 

Only the s-wave contribution to a is relevant in discussing the gamma ray excess coming 
from dark matter annihilation since the velocity of the dark matter in the Galactic Center 
region is relatively low. An exception to this is within the sphere of influence of the Milky 
Way supermassive black hole, but this region corresponds to only a fraction of an arc second 
and is below FGST accuracy. As we see later, a\ is mostly singlet for benchmark points 
A and B. Therefore the s-wave contribution to dark matter annihilation to fermions comes 
mostly from the s-channel diagram involving an exchange of the lightest pseudoscalar a\ 
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given in Fig. [TJ It is approximately given by 



N c k 2 U? f m 2 f ml, / ml 

s lf f Xl ■ /1 f (3.3) 



m 2 



Al 



4lT V 2 (4771^—771^) 5 

where 7V C is the number of fermion colors, U\$ is the (1,/) element of the pseudoscalar 
diagonalizing matrix and m ai is the mass of the lightest pseudoscalar. The s-wave contribu- 
tions from heavier pseudoscalars are suppressed by larger masses as well as smaller mixings 
with the singlet. Moreover, s-channel scalar exchange diagrams are s-wave suppressed, i.e. 
a (xiXi -> hi -»■ //) = 0. 

For benchmark point A, the dark matter mass is m Xl = 7.4 GeV. The physical dark 
matter can be expressed in terms of gauge eigenstates as: 

Xi = 0.0017 B° - 0.0031 W° - 0.0141 h u - 0.0046 h d - 0.0001 h - 0.0008 h t + 0.9999 s. 

We need a light pseudoscalar, O(10) GeV, to get a sizeable annihilation cross section. This 
requires 1% tuning in the parameter space in addition to the tuning needed to make the 
singlino the LSP. The lightest pseudoscalar in the benchmark point is mostly singlet with a 
mixing with other types of pseudoscalar given by 

a x = -0.000002 a u - 0.002193 a d - 0.001203 a - 0.003679 a e + 0.999990 a s , 

with its mass is m ai = 18.7 GeV. 

Having the masses and mixing, we can calculate the total annihilation cross section into 
fermion pairs which gives 

(av) = 4.0 x 10" 26 cm 3 /s (3.4) 

where the hadronic final states cross section is 23% of the total cross section and r pairs 
final state makes up the rest. For benchmark point B given in Table [Til the mass of dark 
matter is m Xl = 7.4 GeV and (av) = 3.0 x 10~ 26 cm 3 /s, with the hadronic final states make 
up 23% of it. The annihillation cross sections given above are within the range of suggested 
cross section for explaining the gamma ray excess in the Galactic Center region given in Ref. 

In this model, dark matter annihilation into SM fermions given in Fig. [T]is also responsible 
for giving the dark matter the correct thermal relic abundance. To show this, we calculate 



the relic abundance which is given by 16 1 

tt Xl h 2 w 2.82 x 10 8 Y^m^/GeV), (3.5) 
10 



where 

Y~ l = 0.264 ^gZm P m Xl {a/x f + 3(6 - \a)/x)} . (3.6) 
In the equation above, mp is the Planck mass and g* is the number of relativistic degrees 



of freedom at freeze-out. The freeze-out epoc 
Tf by Xf = m Xl /Tf, and Xf is determined by 



i Xf is related to the freeze-out temperature 



16| 



Xf = In [0.0764 mp(a + 6b/xf)c(2 + c)m Xl /y/g*Xf\ . (3.7) 

The value of c is usually taken as c = \. Approximating g* to be a ladder function, we 
get that, for both of our benchmark points, the freeze-out epoch is Xf — 21 and the relic 
abundance is 

tt Xl h 2 « 0.1, (3.8) 



which agrees with the cosmologically measured abundance [17( • Since the freeze-out temper- 
ature happens to be around the QCD phase transition temperature, g* varies significantly 
over the change of temperature 18| and the result (13. 8ft can change up to 0(1). However the 
relic density is in the correct ballpark, therefore we do not expect that the correction will 
invalidate our result. An adjustment of parameters can be done when taking into account 
of the variation of g* to get the correct density and annihilation cross section. 

The benchmark points A and B serve as examples to show that in principle this model 
can explain the gamma ray excess in the Galactic Center region. However, the excess could 
also be obtained by some other regions in the parameter space as shown in the Appendix 
IB1 One could do a scan on the parameter space to find the favored region of the model. 

Note that in our relic density calculation, we have neglected possible chargino and 
sfermion contributions coming from resonance and coannihilation effects. This is because 
the charginos have masses (9(100) GeV for all of our benchmark points, and we assume that 
the sfermion masses are at least (9(100) GeV, which is consistent with current LEP bounds. 



IV. DIRECT DETECTION 



Having shown that this model can account for the gamma ray excess in the Galactic 
Center region, we now discuss direct detection of dark matter of this model. In this section, 
we will consider constraints from the search for spin independent, elastic scattering of dark 
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matter off target nuclei. The most relevant contribution for the cross section is given by the 
t-channel scalar exchange diagram with the effective Lagrangian: 

i 

In our benchmark points, the only relevant contribution to dark matter detection comes 
from the lightest scalar and a q can be approximated by 

aq « K ' m « Vl * , (4.2) 
V2v q m 2 hl ' 

where m q is the mass of quark q, v q is the scalar vev associated with quark flavor q, V\ q is the 
(1, q) element of the scalar diagonalizing matrix, and is the mass of the lightest scalar. 



Given the partonic interaction between dark matter and quarks, we can follow Ref. 19] to 
get the effective interaction with nucleons: 

£ e // = f P XiXi PP + fnXiXinn, (4.3) 



where f p and f n are related to a q through the relation [19J 



h,n _ \^ JTq ,{p,n) «£_ ^ 



m q 27 9 ^— ' m 

q=u,d,s H q=c,b,t H 



and (n|m g gg|n) = m n f^ q . Numerically, the fj?' 7 ^ are given by [2 

= 0.020 ± 0.004, f% = 0.026 ± 0.005, f* a = 0.118 ± 0.062 
f£ u = 0.014 ± 0.0043, ft A = 0.036 ± 0.008, ft a = 0.118 ± 0.062, 

while f^g^ is defined by 



(4.5) 



f M = 1 _ J2 /<*»> . (4.6) 

We can approximate / p »s / n since fc s is larger than other /t 9 's and For the purpose 
of comparing the predicted cross section with existing bounds, we evaluate the cross section 
for scattering off a single nucleon. The result can be approximated as 



477^/2 
7T 



a S i « (4.7) 



where m r is nucleon-dark matter reduced mass l/m r = l/m n + l/m xi . 
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We are now ready to show that benchmark point A can explain signals reported by 



CoGeNT {9] and DAM A lOj. For this benchmark point, the lightest scalar mass is 



11.3 GeV. This lightest scalar is mostly singlet and its mixing with other scalars is given by 

h x = 0.089 h u + 0.004 h d + 0.010 h + 0.004 h e + 0.996 h s . 

As in the case of pseudoscalar, contributions from higher mass scalars are suppressed by 
their masses and their mixings with the singlet. The spin independent cross section for the 
benchmark point now can be calculated and is given by 

a S i = 1.7 x 10^ 40 cm 2 , (4.8) 



which is inside the CoGeNT and DAMA favored region [111 ]. 

Similarly, we can show that benchmark point B given in Table [Til has the lightest scalar 
mass rrif ll = 41.5 GeV and spin independent cross section a si = 1.2 x 10~ 42 cm 2 . This 
cross section is two orders of magnitude lower than the present CDMS and XENON bound 
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V. BOUNDS ON THE MODEL 

In this section we discuss various collider bounds that apply to the model. We will spend 
most of the discussions in this section for the benchmark point A given in Table [TJ The 
bounds for benchmark point B as well as the summary of the bounds for benchmark point 
A are given in Table IIII1 

In this model, the decays Z — > X1X1 an d Z —> h\a\ are allowed kinematically. The Z 



decay width has been measured precisely and is given by V = 2.4952 ± 0.0023 GeV [21]. 

Corrections to the decay width can be used as a bound on the mixing between the singlet 
and the Higgs sector. The partial decay width of Z — > X1X1 is given by 

G F 9l . ( 4m 2 \ 2 

r — iwfr-H 1 - it ) ' (5 - 1) 

where Gp is the Fermi constant, rriz is Z mass, and 8 X is given by 

9 X = \W ul \ 2 - \W dl \ 2 + l^oil 2 - \Wa\ 2 ■ (5.2) 
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Benchmark point 


A 


B 


m Xl (GeV) 


7.4 


7.4 


m + (GeV) 


118 


118 


m hl (GeV) 


11.3 


41.5 


m B i (GeV) 


18.7 


19.3 


r^^y.y, (GeV) 


1.4 x 1(T 9 


1.4 x 10" 9 


T z ^h iai (GeV) 


1.1 x KT 11 


4.9 x 10" 12 


it 


8.0 x 1(T 3 


1.3 x 10~ 2 


^modeiCe+e - -thiai) 


1 x 1(T 10 


i x i(r 10 


S mo dei{e + e~ ->■ /i 2 ai) 


i x i(r 12 


2 x 1(T 12 


cr e+e -_^ XlX2 (pb) 


1 x 1(T 5 


1 x 10" 5 



TABLE III: Mass spectrum and bounds for benchmark points A and B. The variable k is given by 
k = (Thz/o'hZ ancl Smodei = ^h iaj I ' &ref , where (Jhidj is the hidj production cross section and a re f 



is the reference cross section defined in Ref. 



2J 



In the equation above, Wfi is the (/, 1) element of the neutralino diagonalizing matrix. The 
decay width of Z — > h\ai is given by 



r 



Z— >h\a\ 



GF\0ha\ 2 3 

-v , 



where 



and 



3\/27r 

6 ha = UuiVui — UdiVdi + f/oiVoi — UeiVei, 

[(m| - (m hl + m ai f) (m| - (m hl - m ai ) 2 )} 



1 



4m| 



(5.3) 



(5.4) 



(5.5) 



For the benchmark point, the partial decay widths in both cases are given by 

T z ^ XlXl = 1.4 x 10~ 9 GeV, 



1.1 x KT 11 GeV 



(5.6) 



-L Z— >h\a\ 

which is well within the measurement error. 

Another bound on the model comes from scalar and pseudoscalar direct production at 
LEP. At LEP a light scalar can be produced by Higgsstrahlung process e + e~ — > Z — > Zh\. 
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Ref. 22[ gives a bound on the coupling strength of Z pairs to scalars regardless of the 



scalar's decay mode. The bound is given in terms of the quantity 



k(m h ) = -U- (5.7) 



SM ■ 
°hZ 



In our model, k(m,h) is given by 



1 



k(m h J = — \v u V ui + v d V di + v V 0i + v e Vu\ , (5 A 

v 

ew 



and its value for the lightest scalar at our benchmark point is 

k(m hl ) = 8.0 x 10~ 3 . (5.9) 
The bound on k(rrih) for the benchmark point hi mass is given by 

fc(11.3 GeV) < 0.09. (5.10) 

Therefore ^(m^) does not exceed the bound from Higgsstrahlung process in our benchmark 
point. The pseudoscalar can also be produced at LEP by the process e + e~ — > Z — > ha. In the 
benchmark point, both h\a\ and h 2 a\ production are kinematically allowed. LEP bounds on 



scalar and pseudoscalar production for various final states are given in Ref. (23|. The bound 
is given in term of S95 = a max /a re f where a max is the largest cross section compatible with 
data and a re f is the standard model hZ production cross section multiplied by a kinematic 
scaling factor. Defining S mo d e i = o'hm/ /(7 refj where <Jh iaj is the model's /ijOj production cross 
section, the bound on the model is given by S mo( iei < Sg 5 . For our benchmark point, S mo d e i 
is given by 

Smodei (e + e~ Mi) = 1 x 10~ 10 , 

(5.11) 

Smodei (e + e -> h 2 a 1 ) = 1 x 10 , 

which is lower than the bound, S95 ~ O(10~ 2 ), in both cases. 

We note that the lightest chargino mass is 118 GeV for the benchmark point, which 



exceeds the PDG bound of 94 GeV 21|. In the case of a long lived chargino however, the 
bound can be made much stronger and is currently at 171 GeV. We have calculated the 
lifetime of the chargino in our model assuming a stau mass of 110 GeV and have found that 
it is short lived, thus this latter bound is not of concern. We should point out however, that 
our analysis has been done at tree level. Loop corrections could change these results but are 
beyond the scope of this paper. 

15 



Finally, we need to calculate the bound on neutralino productions. Ref. [24J discusses the 
bound on production of the lightest and second to lightest neutralinos at LEP, e + e~ — > X1X2, 
where \2 decays into Xiff- Assuming that the selectron is much heavier than the Z, the 
main contribution comes from s-channel Z exchange. For our benchmark point, we calculate 
the cross section to be 

a e+e -^ xlX2 = 1 x 10~ 5 pb, (5.12) 

while the bound is O(0. 1) pb. A summary of all these bounds is given in Table HTT1 

The light particles are mostly singlet and have very little mixing with the Higgs sector. 
This make the particles unlikely to be produced at near future experiments. However the 
heavier sector has a richer phenomenology. For example, heavier scalars are mostly h u , hj, 
ho, and he therefore they have a better chance of being detected in future colliders {§]. 



VI. CONCLUSIONS 



In this paper, we have presented a supersymmetric model of 7 — 10 GeV dark matter, 
which is capable of describing the FGST observations. In a recent analysis of FGST data, 
Hooper and Goodenough found an excess in gamma ray emission from within 1.25° of the 
Galactic Center. They showed that this can be explained by annihilating dark matter if 
the dark matter has a mass between 7—10 GeV, annihilates into r-pairs most of the time, 
but into hadronic channels the other 15 — 40% of the time, and (av) falls within the range 
4.6 x 10~ 27 — 5.3 x 10~ 26 cm 3 /s [l|. Our model achieves these requirements by minimally 
extending the SLHM to include a scalar singlet whose superpartner is the dark matter 
particle. Due to the Yukawa structure of the SLHM the scalar particles mediating the dark 
matter annihilation have an enhanced coupling to leptons. This provides a natural means 
for satisfying the second requirement put forward by Hooper and Goodenough. 

We have shown that this model produces the correct dark matter thermal relic density 
and is consistent with current collider bounds. In addition, we have shown that this model 
is consistent with the direct detection signals reported by both CoGeNT and DAMA for 
certain regions of parameter space, while for other regions of parameter space, the model 
yields a spin independent cross section far below the present CDMS bound, but maintains 
the right relic density and continues to explain the FGST observations. Thus our model is 
fully able to accommodate the results reported by CoGeNT and DAMA in the case of their 
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Field 


^3a 


^•3£ 


Field 


^3o 




H u 


UJ 


1 


-^01 


1 


1 


H d 


OJ 


1 


Xq2 


a; 2 


a; 2 


H 


1 


U) 


X q l 




1 


H e 


1 


UJ 


X q 2 


w 2 


1 


E 


1 


UJ 2 


Xg\ 


1 


u> 


Q 




1 


Xg2 


1 


UJ 2 



TABLE IV: Transformation rule for the Z^ q x TL^g symmetry. Each field transforms as (f> — > X(f>, 
where X is the corresponding factor shown in the table. For each case, u) = 1. Other fields not 
shown in the table are neutral under Z^ q x /L^g 

vindication, but it is in no way contingent upon their validity. 
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Appendix A: Breaking Terms 

In this appendix, we discuss a possible source of the terms in V^ ft that break the Z 2 
symmetry of the superpotential. Generally, one can imagine such breaking terms arising 
from the F-term of some hidden sector superfield receiving a vacuum expectation value. To 
be more specific, we consider a possible scenario that results in such breaking terms and also 
explains the smallness of K q and Kg. In this scenario there is a hidden sector, which contains 
the six fields X i, X 02 , X q i, X q2 , Xg\ and Xg 2 . The F-terms of the fields receive vevs 

(F^) ~ O(10 n GeV) 2 , (Al) 

so that 

M SUSY ~ (A2) 
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is at the TeV scale. The index i denotes 01, 02, ql, q2, £1, and £2. A Z^ q x Z 3 £ symmetry 
is imposed, under which the fields transform according to Table IIVI The hidden sector 
fields Xi couple to visible sector fields in a high energy, fundamental theory, and are Planck 
suppressed in the low energy effective theory. Consequentially, the lagrangian contains terms 
such as 

AC = W S d ^iXo2H u H £ + |L J d 2 6X 02 SH u H e + h.c, (A3) 

where d 2 9 = d{99) and d A 6 = d(66)d(60) represent integration over Grassmann variables 
and /' and m! are coupling constants. When the F-terms of Xq\ and Xq 2 receive vevs, the 
terms in Eq. (1A3|) give rise to 



AC = f{F °^ p F02) Jd%mmH u H e + ^^- Jd 2 9(99)SH u H £ + h.c. 



f(Foi)(F 02 ) u u , m'(F 02 ) 



H u H e + SH u H e + h.c. (A4) 



Ml ' M P 



-> nlH u H e + fi c SH u H e + h.c. 
Similarly, the breaking parameters [x\ and fi d arise from the Planck suppressed terms 

AC= M 2 ~ J d49 ^AH H d +^f p J d 2 6X 02 SH H d + h.c. 



a' 



^ (F °;]f 02) H H d + SH H d + h.c. (A5) 

Mi, Mc 



lp lVlp 

.2 ; 



-»■ fJ. 4 H H d + fi d SH H d + h.c, 
while the parameters m^ and arise from 

ac = jp J d 4 ex^x n H^H + Jl y d^ex^XnHim + h.c. 

. h'(F Q2 )(Fg) t i / (F 02 )(F a ) t (A6) 

-> rn 2 u0 HlH + m^i^ + h.c. 

In this way, all of the Z 2 breaking terms are generated. At this point it should be noted 
that the Z 3(? x Z 3 ^ symmetry actually prohibits the terms fi q H u H d , figHoHi, K q SH u H d , and 
kiSHqHi from appearing in the superpotential [see Eq. (12.11) ]. As far as the ji q and fie terms 
are concerned, this is not a problem since they are generated by the vevs of the X q2 and Xp 2 
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M P 

V 
Alp 



b 'r- f d^OX^HvHi + h.c. 



£jd*exj 1 sp + h.c. 

\r J d 4 (d'X \x ql + <t'x£x 01 + d'"X \x e2 + d""X q \x q2 ) H u H d + h.c. 
\r J d 4 9 (e'X [x a + e"X e lx 01 + e'"X \x q2 + e""X e \x e2 ) H H e + h.c. 



M 
Jl 



j^Jd*6 (f'X \x 02 + f"X q \x e2 + f"X e {x q2 ) H u H e + h.c. 



jjjr / d 4 ^ (3^01^02 + + gPX&Xgz) H H d + h.c. 

/ d 4 6 (h'X&Xti + ^"^1^02 + h"'X q \x i2 + /i""^X ffl ) fl-Jflo + h.c. 
\r J d*6 (tX&Xti + i"X q \x Q 2 + i"'X^X a + i""X £ [x ql ) Hj H e + h.c. 



M 
M 



Al 



kJd 4 9^ i j i X i 1 X i H}H f + h.c. 
£ f d 2 9X ql SH u H d + h.c. 



M 



-fefdPOXtiSHoHi + h.c. 
jfe J d 2 eX 02 SH u H e + h.c. 
^ / d 2 9X 02 SH H d + h.c. 

JT/^^^X^ + h.C. 



5 g ? /d 2 0X o 5 3 + h.c. 



/ d 2 9fi q H u H d + h.c. 
/ d 2 6^ i H H i + h.c. 
Jd 2 6X 2 S 2 + h.c. 
^?flufr d + h.c. 
H 2 HqHz + h.c. 
H 2 H u Hi + h.c. 
H 2 H H d + h.c. 
m 2 u0 HlH + h.c. 

+ h -c- 
mj|-ff/| 2 + h.c. 
/j, a SH u H d + h.c. 
HbSHoHg + h.c. 
H c SH u H e + h.c. 
H d SH H d + h.c. 

&!?S 2 +h.c. 

a s 5 3 + h.c. 



TABLE V: A complete list of superpotential and V^ Q ft terms generated by the Xi in this example. 



fields in the same manner 

a' 



AC 



a'(F 32 ) 



d 4 9X q2 H u H d 



V 



d 2 6d 2 9(99)H u H d 



d 4 9Xg2 Hq Hi 
M P 



d 2 9d 2 9(99)H H l 



d 2 9H u H d 



b'{F e2 ) 



(A7) 



c( 2 9HqHi 



M P J M P 
^fi q J d 2 9H u H d + ^ J d 2 9H H e . 

In this UV completion scenario, the terms corresponding to « 9 , Ke, Ai and i are not generated 
in this way. Because of the Z 3g x Z 3 ^ symmetry, they are entirely absent at tree level. 
Benchmark points II and V in Table [VI] satisfy n q = Ki = Ai = t = and yield results 
consistent with our goals. Since we are not committing to this particular UV completion 
scheme, we consider several other benchmark points that include nonzero values for these 
parameters. A list of the soft breaking terms relevant to this paper, which are generated by 
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the fields X iy is given in Table IVl 

Appendix B: List of benchmark points 

In this Appendix, we show several benchmark points given in Table IVIl Benchmarks 
point I-III lie in the suggested CoGeNT and DAMA range, while benchmarks point IV- VI 
satisfy CDMS bound. Benchmark point I is identical with benchmark point A discussed in 
the text. Benchmark point IV is identical with benchmark point B. Benchmark points II 
and V are motivated by mechanism described in Appendix [A] 
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TABLE VI: Additional benchmark points 



Benchmark point 


I 


II 


III 


IV 


V 


VI 




0.01 





0.01 


0.01 





0.01 


i 


0.01 





0.01 


0.01 





0.01 




0.6 


0.6 


0.5 


0.6 


0.6 


0.5 


tan a 


20 


15 


30 


20 


30 


25 


tan/3 


50 


30 


30 


50 


25 


25 


tan 


10 


10 


5 


10 


5 


5 


v s (GeV) 


50 


50 


100 


50 


50 


100 


v u (GeV) 


245.6 


245.3 


245.7 


245.6 


245.7 


245.6 


v d (GeV) 


4.9 


8.2 


8.2 


4.9 


9.8 


9.8 


v (GeV) 


12.2 


16.2 


8.0 


12.2 


8.0 


9.6 


v t (GeV) 


1.2 


1.6 


1.6 


1.2 


1.6 


1.9 


u (GeV) 


125 


125 


200 


125 


125 


150 


M/ (GeV) 


125 


125 


150 


125 


150 


150 


A? (GeV 2 ) 


100 2 





150 2 


100 2 





50 2 


A 9 (GeV) 


-35 


-35 


-63 


-35 


-35 


-63 


Mi (GeV) 


500 


500 


250 


500 


250 


200 


Mo (GeV) 


500 


500 


500 


500 


500 


400 


m 2 n (GeV 2 ) 


-100 2 


-150 2 


-150 2 


-100 2 


-150 2 


-150 2 


m 2 (GeV 2 ) 


100 2 


200 2 


100 2 


100 2 


200 2 


100 2 


u? (GeV 2 ) 


400 2 


300 2 


300 2 


400 2 


400 2 


350 2 


M9 (GeV 2 ) 


200 2 


300 2 


250 2 


200 2 


200 2 


300 2 


Ml (GeV 2 ) 


200 2 


200 2 


250 2 


200 2 


250 2 


200 2 


u 2 a (GeV 2 ) 


400 2 


200 2 


200 2 


400 2 


400 2 


100 2 


Ma (GeV) 


100 


75 


75 


100 


100 


80 


M6 (GeV) 


200 


150 


300 


200 


250 


400 


/i c (GeV) 


200 


200 


400 


200 


300 


200 


Md (GeV) 


200 


100 


100 


200 


250 


100 


Continued on the next 


page 
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TABLE VI: continued 



Benchmark point 


I 


II 


III 


IV 


V 


VT 
V 1 


t (Lrev J 


60. 6 3 









83.9 3 




55.0 3 









— o7.y 


b s (GeV ) 


63. 4 2 




43. 6 2 




98. 2 2 




66. 3 2 




47.1 2 




on n2 

yy.u 


a s (GeV) 


—49 4 




91 7 




9 




—49 9 




9n n 




—50.2 


m Xl (GeV) 


7.4 




7.4 




7.7 




7.4 




7.4 




7.7 


m ± (GeV) 


118 




117 




151 




118 




117 




137 


m hl (GeV) 


11.3 




19.2 




12.8 




41.5 




41.4 




26.1 


m ai (UeV J 


18.7 




16.1 




18.8 




19.3 




19.2 




ii 7 
11. / 


/_„,\ /cm 3 \ 


4.0 x 10" 


-26 


3.4 x 10" 


-26 


4.6 x 10" 


-26 


3.0 x 10" 


-26 


3.1 x 10" 


-26 


A i \y 1 n— 26 
4.1 X 1U 


{crv (xiXi~ >hadrons)) 
(av) 


23% 




38% 




32% 




23% 




24% 




30% 


o"5/( cm z ) 


1.7 x 10" 


-40 


1.2x10" 


40 


1.5 x 10" 


-40 


1.2 x 10" 


-42 


6.1 x 10" 


-42 


1.5 x 10 -41 


T Z^ X iXi ( GeV ) 


1.4 x 10 


-9 


n 




2.1 x 10" 


-10 


1.4 x 10 


-9 


n 




c o \ , in — 10 

6.3 x 10 


rz^/! iai (GeV) 


1.1 x 10" 


-11 


1.2 x 10" 


-10 


1.4 x 10" 


-10 


4.9 x 10" 


-12 


4.2 x 10" 


-11 


1.2 x 10" 10 


k 


8.0 x 10 


-3 


3.5 x 10 


-2 


2.2 x 10 


-2 


1.3 x 10 


-2 


0.12 




2.8 x 10" 2 


S mo dei(e + e~ -)• hid) 


1 x 10" 


10 


2 x 10" 


9 


2 x 10" 


9 


1 x 10" 


10 


1 x 10" 


9 


2 x 10" 9 


5'mod e ;(e + e~ ->■ /l20l) 


1 x 10" 


12 


5 x 10" 


11 


3 x 10" 


11 


2 x 10" 


12 


1 x 10" 


10 


4 x 10" 11 


cr e+e -^. XlX2 (pb) 


1 x 10" 


5 







5 x 10" 


9 


1 x 10" 


5 







4 x 10" 6 
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